Coupling of phonons and electromagnons in GdMnOs 
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The infrared and Terahertz properties of GdMnOs have been investigated as function of tempera- 
ture and magnetic field, with special emphasis on the phase boundary between the incommensurate 
and the canted antiferromagnetic structures. The heterogeneous incommensurate phase reveals 
strong magnetodielectric effects, characterized by significant magnetoelectric contributions to the 
static dielectric permittivity and by the existence of electrically excited magnons (electromagnons) . 
In the commensurate canted antiferromagnetic phase the magnetoelectric contributions to the dielec- 
tric constant and electromagnons are suppressed. The corresponding spectral weight is transferred 
to the lowest lattice vibration demonstrating the strong coupling of phonons with electromagnons. 



Multiferroic materials with the simultaneous occur- 
rence of magnetism and ferroelectricity, are a hot topic in 
recent solid-state research. They provide interesting and 
spectacular physical properties and promise attractive 
applications P, ■ Multiferroic behavior occurs in a 
variety of systems originating from very different physical 
mechanisms, including materials with independent mag- 
netic and ferroelectric subsystems, like some boracites, 
Aurivillius phases, hexagonal manganites, and the lone- 
pair ferroelectrics with magnetic ions Q]- Recently a 
new class of multiferroics, namely ferromagnetic sulfos- 
pinels with relaxor-like short range ferroelectric (FE) or- 
der have been detected 3 with a strong coupling of the 
electric and magnetic properties al low frequencies. In 
these spinel compounds the ferromagnetism is induced 
via strong indirect exchange interaction, but the origin 
of ferroelectricity remains unclear sofar. Finally, in the 
perovskite manganites there is robust experimental evi- 
dence 0, 13 that the onset of helical magnetic order in- 
duces spontaneous FE polarization 0,13 ■ Dzyaloshinskii- 
Moriya type interactions have been utilized to explain the 
ferroelectricity which is induced by the helical spin struc- 
ture llOl lllj l . A similar spin-driven ferroelectricity is 
believed to be operative in Ni3V208 12]. 

After having established the ground-state proper- 
ties of this interesting class of materials, the study of 
their novel dynamic properties will significantly enhance 
our knowledge of the magneto-electric (ME) coupling 
. Magnons are the characteristic excitations of mag- 
netic structures, while soft phonons as inferred by the 
Lyddane-Sachs- Teller relation condense at canonical fer- 
roelectric phase transitions. It seems clear that soft 
phonons cannot be relevant excitations in the ferroelec- 
tric manganites, as (improper) ferroelectricity is induced 
by the magnetic order coupled to the lattice. Recently it 
has been shown that electro-magnons are relevant collec- 
tive modes in this new class of ferroelectrics . Electro- 
magnons are spin waves that are excited by an ac electric 
field. By measurements in TbMnOa and GdMnOa it has 
been documented that these new excitations exist not 
only in the magnetic phase characterized by the helical 
spin structure, but also in the longitudinally modulated 



(sinusoidal) structure, provided that a " helical- type" vec- 
tor component of the spin- wave is dynamically induced 
via the ac electric field jlJI . 

In this Letter we present detailed investigations of the 
Terahertz and FIR properties of GdMnOs. We investi- 
gate electromagnons and phonons as function of temper- 
ature and magnetic field. We provide striking experimen- 
tal evidence that i) electromagnons are strongly coupled 
to phonons with a considerable shift of optical weight be- 
tween these excitations and ii) electromagnons contribute 
considerably to the static dielectric constant. 

Single crystals of GdMnOs have been prepared using 
the floating-zone method with radiation heating. The 
samples have been characterized using X-ray, magnetic 
and dielectric measurement |l5j[ . The basic properties of 
our sample s ag ree well with the results obtained by other 
groups la, Il6| . The experiments at Terahertz frequen- 
cies were carried out in a Mach-Zehnder interferometer 
|l7j which allows the measurements of transmittance and 
phase shift. The absolute values of the complex dielectric 
permittivity e* = ei+ is2 were determined directly from 
the measured spectra using the Fresnel optical formulas 
for the complex transmission coefficient. The spectra in 
the infrared frequency range have been obtained using 
a Bruker IFS-113 Fourier-transform spectrometer. The 
experiments in external magnetic fields were performed 
in a superconducting split-coil magnet with polyethylene 
windows allowing to carry out reflectance experiments in 
magnetic fields up to 7T. 

The upper panel of Fig. ^ shows the H-T phase dia- 
gram of GdMnOs for H ||c in the zero-field-cooled (ZFC) 
regime ^3 which basically coincide with the diagrams 
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published previously 0, 0| . GdMnOs is paramagnetic 
above Tn — 42 K and transforms into an incommensu- 
rate antiferromagnetic (IC-AFM) state below this tem- 
perature. Depending upon the value of the external 
magnetic field along the c-axis, the magnetic state be- 
comes canted antiferromagnetic (CA-AFM) below 20 K. 
For r < 9 K and in low fields the spin structure re- 
veals increasing complexity due to an additional order- 
ing of Gd-subsystem. This region, which is indicated by 
a hatched area in the phase diagram of the upper panel 
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FIG. 1: (color online) Upper panel; Magnetic phase diagram 
of GdMnOa for H\\c for zero-field cooling (ZFC). Circles and 
stars: results of magnetization measurements, triangles: data 
from dielectric experiments [3 [ . Lines are drawn to guide the 
eye. PM - paramagnetic, IC-AFM - incommensurate antifer- 
romagnetic (probably sinusoidally modulated 0), CA-AFM 
- canted antiferromagnetic structure. The hatched region in- 
dicates the ordering of Gd-sublattice. 

Lower panel: Reflectance of GdMnOa with polarization along 
the a-axis and at phonon frequencies. Symbols: experi- 
ment, solid line: model calculation assuming the sum of 15 
Lorentzians (see text). 



of Fig. n is not further discussed in the course of this 
work. There is one important difference compared to the 
{H, T) phase diagram published in Ref. (15| : under zero- 
field cooling conditions no phase transition occurs and 
the AC-AFM state remains stable down to 9 K. This 
fact allows to switch between the magnetoelectric IC- 
AFM and the CA-AFM phases with very low fields. The 
IC-AFM region is especially interesting from the spec- 
troscopic point of view, because unusual excitations of 
mixed magnetoelectric origin (electromagnons) exist . 
As have been shown previously, the electromagnons are 
magnons which can be excited by the electric component 
of the electromagnetic wave. These excitations are sup- 
pressed in the CA-AFM state. It is the aim of this work 
to study magnetic-field and temperature dependence of 
electromagnons in a broad frequency range and to inves- 
tigate their possible coupling to phonon modes. 
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FIG. 2: (color online) Magnetic field dependence of the pa- 
rameters of the phonon at 119 cm~^ crossing the boundary 
between IC-AFM and CA-AFM phases. Upper panel: dielec- 
tric contribution (Ae), lower panel: eigenfrequency (uq). The 
lower inset shows the reflectance spectra at typical magnetic 
fields. The upper inset shows the changes in dielectric permit- 
tivity at a frequency close to electromagnon, which demon- 
strates direct connection of both effects. 



The lower panel of Fig. ^ shows the reflectance spec- 
trum of GdMnOs at phonon frequencies and for the ac 
electric field component e parallel to the crystallographic 
a-axis. We note that this direction of the electric field 
reveals large effects in the temperature and field de- 
pendence of the dielectric constant and is especially 
important for the magnetodielectric effects in this com- 
pound. The solid line in the lower panel of Fig. ^ has 
been calculated using a sum of 15 Lorentzians. Here 
ten strongest excitations correspond to phonons polar- 
ized along the a-axis at frequencies 119, 188, 231, 308, 
325, 400, 460, 475, 509, and 568 cm"!. Weaker features 
at 296, 441, 539, and 641 cm^^ are due to leakage of 
the b-axis component, most probably because of a small 
misorientation of the crystal. The weak excitation at 76 
cm~^ is of unknown origin and might be due to a forbid- 
den exchange excitations of Mn or Gd subsystem, which 
becomes IR allowed via strong magnetoelastic coupling. 

Except for the lowest frequency phonon at 119 cm~^, 
no measurable changes in the phonon parameters have 
been detected between different magnetic phases. On 
the contrary, the 119 cm~^ phonon, which is observed 
for a-axis polarization only, reveals substantial changes 
between the IC-AFM and CA-AFM phases. We recall 
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that the peculiarity of this transition into the CA-AFM 
phase is the occurrence of strong magnetodielectric ef- 
fects, e.g. the magnetic field dependence of the a-axis 
dielectric constant in the frequency range from zero to 
about one Terahertz (~ 40 cm~^) As have been 

shown previously 14], these effects can be directly re- 
lated to the existence of electromagnons at ~ 20 cm~^. 

The magnetic field dependence of the 119 cm~^ phonon 
excitation on crossing the IC-CA magnetic phase bound- 
ary is shown in the lower inset of Fig. |21 which repre- 
sents the reflectance spectra at three typical fields. The 
dielectric contribution (Ae, strength of the mode) and 
the eigenfrequency (vq) of the phonon are shown in the 
main panels. The transition is manifested by a signifi- 
cant softening of the phonon frequency and by a 60% in- 
crease of the mode strength on increasing magnetic field. 
The characteristic fields and the width of the transition 
closely coincide with the corresponding changes of the 
dielectric constant at lower frequency, which is demon- 
strated in the upper inset of Fig. El The difference in ei 
on increasing and decreasing fields documents the meta- 
stable character of the IC-CA phase boundary. Similar 
behavior of the high- and low-frequency dielectric con- 
stant indicates that the same mechanism is operative for 
both effects, namely that at the IC-CA transition elec- 
tromagnons are suppressed by the external magnetic field 

In order to analyze the interplay between electro- 
magnons and phonons the complex dielectric permittiv- 
ity has been calculated from the reflectance spectra via 
the Kramers-Kronig transformation adding the Terahertz 
spectra at low frequencies. The results both in IC-AFM 
(15 K, T) and in CA-AFM state (15 K, 2 T) are shown 
in Fig. O Here the data above 40 cm~^ represent the re- 
sults of the Kramers-Kronig analysis and the data below 
this frequency have been taken directly from the Ter- 
ahertz transmittance experiments. The lower panel of 
Fig. O clearly demonstrates the over-damped almost re- 
laxational nature of the electromagnon and its suppres- 
sion by the external magnetic field. On the other hand, 
we know from Fig. |21that the phonon mode at 119 cm~^ 
gains considerable spectral weight on increasing magnetic 
field. The substantial spectral weight which is removed 
from the low-frequency range, with the dielectric con- 
stant decreased for all frequencies below ~ 40 cm^^, is 
transformed into phonon intensity at 119 cm~^. In order 
to obtain the estimate of the spectral-weight transfer, the 
complex dielectric permittivity has been fitted assuming 
the phonon parameters as obtained from the fits of the re- 
flectance and a single over-damped Lorentzian represent- 
ing the electromagnon with — 40 cm^^, 7 = 75 cm^^ 
and Ae — 1.6. Here 7 is the damping. These parameters 
correspond to the maximum in 62 at ^'max — 27 cm~^ 
and to the spectral weight S = Aevq ~ 2.5 • lO'^ cm^^. 
From the data in Fig. [21 the spectral-weight increase of 
the phonon is obtained as AS{IC-CA) ~ 1.7-10^ cm"^. 
Taking into account strong scattering of the data between 
40 and 100 cm~^ we conclude rough coincidence of these 
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FIG. 3: (color online) Spectra of the a-axis dielectric per- 
mittivity of GdMnOa in the frequency range between electro- 
magnon at 25 cm~^ and phonon at 119 cm~^ at T and 2 
T. Upper panel: real part of the dielectric permittivity, lower 
panel: imaginary part. Triangles represent the results from 
the Terahertz transmittance. Circles - spectra obtained via 
Kramers-Kronig analysis of the reflectance. Solid line corre- 
sponds to the sum of Lorentzians as described in the text. 

values and the transformation of the spectral weight of 
the electromagnon into the phonon mode at 119 cm~^. 

A remaining question is why electromagnons are ob- 
served both in GdMnOa, in a non- ferroelectric colhnear 
magnetic structure, as well as in TbMnOa, in a ferroelec- 
tric with helical spin structure jl4j . In order to under- 
stand this point. Fig. 01 provides a closer inspection of 
the free-ME energy in terms of the Fourier components 
of the dynamic variables Aj, related to the main AFM 
vector A of modulated spin structure |18| (see also Refs. 

mil): 

k 

-zP.^ar(A, X A^).^ 

k 

= -^PE„t(Afe) , (1) 

k 

where P is electric polarization. The ME coefficients 
^a;z,2j: ^j^^ nearest neighbors within the ab-plane 
are determined by ti^^'^^ = 2Afa^^^^^ cos(27rfc6) sin(7rfc6), 
where a^^'^^ are constants, N the number of Mn ions. 
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FIG. 4: (color online) Schematic presentation of magnetic 
structures, showing different magneto-electric behavior, (a) 
Homogeneously modulated structure: no static electric po- 
larization (APme ~ 0), no ME contribution to the dielec- 
tric permittivity (Ae — 0). (b) Inhomogeneously modulated 
coUinear structure: no electric polarization, but nonzero con- 
tribution to the dielectric permittivity, (c) Cycloidally mod- 
ulated structure: finite electric polarization and ME contri- 
bution to the dielectric permittivity. 



and b the lattice constant. This expression was de- 
rived using the crystallographic symmetry D^f^ (Pbnm 
space group) and a modulated spin structure with k ~ 
(0, fc, 0). We omitted weak contributions from other AFM 
vectors F, C, G which exist in this structure. In 
space representation and continuum limit Eq. Q cor- 
responds to $„je = -a^Px{AxdAy/dy - AydA^/dy) - 
a,P,iA,dAy/dy ~ AydAJdy) iW], and for = 
is reduced to Dzyaloshinski-Moriya type interactions 

It is clear that in a homogeneous magnetic state, like in 
the CA-AFM phase (Fig. 4a), the ME free energy is zero 
and no contribution to the dielectric constant, no electro- 
magnons and no spontaneous polarization can exist. To 
obtain the ME contribution to the electric susceptibility 
in sinusoidal and spiral states we consider the total free 
energy of the system: 

$(Afc,P) = 



-E • P + PV2xs + - TS{Ku) , 



(2) 



where the first three terms correspond to exchange and 
anisotropy energies, the forth and fifth terms represent 
dielectric contributions from external electric fields E and 
the last term is the spin entropy. By minimizing Eq. (j^J 
with respect to P the free energy can be represented as 
a function of the non-equilibrium values of Ak- In a si- 
nusoidal spin structure with A^ — (0, A^, 0) and keeping 
only the main harmonic fcg of the modulated structure, 
the electric susceptibility, e.g. along the a - axis, can be 
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dm 



92$ 



dE^dAl dE^dAl* ' dA^dAl* 



XE + {xEaro?Al,Ai;/Kta 



(3) 



In the sinusoidal phase no spontaneous polarization can 
exist since only one A^ component is nonzero, but the 
ME contribution to electric susceptibility arises accord- 
ing to Eq. 0. It originates from an electric field- 
induced rotation of the spins in the ab-plane, i.e. from 
A^ spin components. Similar contribution can also ex- 
ist along c-axis. Their relative values are determined 
by both, the corresponding magnetoelectric (a^^'^^) and 
magnetic anisotropy (Kba^tc) constants. Finally in case 
c) of Fig. 0] with helical or cycloidally modulated 
spins, e.g. A^ = {0, Af, A^), a spontaneous ferroelec- 
tric polarization along the c-axis Pz, finite contribu- 
tion to the dielectric constant along the a-axis Xe ~ 
Xe + [XEd^QY S'^ I Kbai and electromagnons exist. 

Very recently Katsura et al. J^J have calculated the 
collective mode dynamics of helical magnets coupled to 
the electric polarization. For the ac dielectric proper- 
ties their main findings are the occurrence of two modes, 
one of which is derived from the phonon mode with a 
frequency close to the eigenfrequcncy of the uncoupled 
phonon, and one originating from the spin wave with 
a frequency proportional to V SJD, where S is the spin 
value, J the exchange coupling and D the anisotropy. Us- 
ing realistic parameters they calculate an electromagnon 
frequency of 10 cm^^, close to the experimental obser- 
vation. It also follows from their calculation that in the 
electromagnetic phase the phonon eigenfrequency (i^o) is 
enhanced by v^j^v^, where is the frequency of the 
spin-plane rotation mode, which from the lower panel of 
Fig. 2 can now be estimated as 13 cm^^. 

In conclusion, studying the low-frequency electrody- 
namics of GdMnOs with a finite magneto-electric (ME) 
coupling, we were able to demonstrate i) the existence 
of electromagnons, ii) that these collective modes of ME 
magnets contribute to the static dielectric constant, and 
iii) that at the transition to a homogeneous magnetic 
phase with no ME coupling the electromagnons are wiped 
out and their spectral weight is transferred to an opti- 
cal phonon, which in addition reveals a slight softening 
of the eigenfrequency. To conclude with an outlook for 
the future applications, a material with similar effects at 
room temperature will directly allow the design of a new 
class of magneto-optical devices. The dielectric constant 
and hence the refractive index n ~ ^j£\ can be tuned by 
moderate magnetic fields. 

This work was supported by BMBF(13N6917B - 
EKM), by DFG (SFB 484), and by RFBR (04-02-16592, 
06-02-17514). 
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